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Network data analysis methods are the only way to properly separate real gravitational wave 
(GW) transient events from detector noise. They can be divided into two generic classes: the coin- 
cidence method and the coherent analysis. The former uses lists of selected events provided by each 
interferometer belonging to the network and tries to correlate them in time to identify a physical 
signal. Instead of this binary treatment of detector outputs (signal present or absent), the latter 
. . . method involves first the merging of the interferometer data and looks for a common pattern, consis- 

, tent with an assumed GW waveform and a given source location in the sky. The thresholds are only 

. applied later, to validate or not the hypothesis made. As coherent algorithms use a more complete 

' information than coincidence methods, they are expected to provide better detection performances, 

' ^ , but at a higher computational cost. An efficient filter must yield a good compromise between a low 

CL(' false alarm rate (hence triggering on data at a manageable rate) and a high detection efficiency. 

D ' Therefore, the comparison of the two approaches is achieved using so-called Receiving Operating 

Characteristics (ROC), giving the relationship between the false alarm rate and the detection effi- 
QQ ' ciency for a given method. This paper investigates this question via Monte-Carlo simulations, using 

I the network model developed in a previous article. Its main conclusions are the following. First, 

. a three-interferometer network like Virgo-LIGO is found to be too small to reach good detection 

^ ' efficiencies at low false alarm rates: larger configurations are suitable to reach a confidence level high 

, enough to validate as true GW a detected event. In addition, an efficient network must contain in- 

■ terferometers with comparable sensitivities: studying the three-interferometer LIGO network shows 

J"* ' that the 2-km interferometer with half sensitivity leads to a strong reduction of performances as 

^— , compared to a network of three interferometers with full sensitivity. Finally, it is shown that coher- 

' ent analyses are feasible for burst searches and are clearly more efficient than coincidence strategies. 

I Therefore, developing such methods should be an important goal of a worldwide collaborative data 

— , analysis. 
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I. INTRODUCTION 



■ The first generation of large interferometric gravitational wave (GW) detectors [1-5] is producing a growing set 
" " of experimental results, showing that the detectors already built are coming close to their foreseen sensitivities. 
Therefore, it is very important to consider exchanging data and perform analysis in common. Such network data 
analysis methods are compulsory to separate - with a sufficient confidence level - real GW transient signals from 
noise occurring in one particular detector. 

GW burst signals have usually a small duration (a few ms) and a poorly known shape: for instance, type II 
supernovae or the merging phase of coalescing compact binary systems belong to this category. As they are not 
accurately modeled, only suboptimal methods [6-11] can be used to detect them. Therefore, various related outputs 
coming from a set of interferometers are required to reach a definite conclusion on the reality of the GW event. 
In addition to the definition of efficient filters suitable to analyze single interferometer outputs, it is also necessary 
to estimate the performances of different network data analysis methods. Studying this problem with Monte-Garlo 
simulations is the goal of this paper. 

Network data analysis methods can be classified into two categories: coincidence or coherent filtering. The first 
method is simpler to use and has been already considered for years. In this approach, each interferometer belonging to 
the network analyzes separately its own data and produces a list of selected events, characterized by their timing and 
their maximum signal-to-noise ratio (SNR) exceeding some suitable threshold levels. In a second step, these events 
are correlated with those found by other detectors in order to see if some are compatible with a real GW source. Many 
articles in the literature [12-14] deal with this topic, in particular with real data taken by resonant bar experiments 
[15,16]. 
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As pointed out independently [17,18], such coincidence analyzes are not optimal in the sense that their binary use 
of interferometer data (a GW signal in a given detector is either present or absent) leaves aside important information 
on the possible correlations between the different datasets. Indeed, as the detector beam patterns are not uniform 
[19], the interaction between the physical signal and a given interferometer depends on the relative location of the 
source. Thus, the GW amplitude scales differently in the various components of the network causing, in addition to 
noise fluctuations, the filtered outputs to exceed or not the selection threshold. These large variations of the GW 
response could suppress even a strong GW signal in a given interferometer [14]. In a coherent analysis, all detector 
outputs contribute to the filtering algorithm. Triggered and merged in a suitable way, a higher statistical significance 
can be achieved compared to individual analyzes, thus improving the network detection potential. 

Coherent methods have already been studied in the literature, particularly for the search of signals with known 
waveforms. Both Ref. [17] and [18] consider network data analysis methods based on a likelihood function, which 
turn out to be direct extensions of the Wiener filtering, optimal for a single detector search when the signal shape 
is known. Coincidence and coherent detections have been compared in a two-detector network [18], with a simple 
model describing its interaction with a GW burst signal (two sinusoid cycles). The coherent search was found to be 
always better than the coincidence method and the results are robust with respect to the noise statistics. Ref. [17], 
later extended in [20] , deals with in-spiral waveforms at Newtonian order. A new formalism is developed based on a 
likelihood function, in a way very similar to matched filtering, but now with a parameter space containing two more 
unknowns which correspond to the source position in the sky: consequently, the data analysis procedure becomes 
more computationally expensive. Ref. [21,22] applied then this framework to post-Newtonian inspiral waveforms and 
to the case of a non-Gaussian noise. 

The present study takes advantage of these pioneering works and applies the same kind of methods for GW bursts, 
using the network detection model originally developed in [14] to study coincidences between interferometers. On 
the one hand, this framework is simple enough to perform a large number of Monte-Carlo simulations, and thus 
to compare accurately coincident and coherent detections. On the other hand, all the features characterizing the 
interaction between a GW and a network of interferometers arc properly taken into account: non uniform angular 
pattern of the detectors, propagation time delays between them and data sampling. 

As an efficient detection algorithm must be a good compromise between a high detection probability and a low 
false alarm rate, a standard tool to estimate the performances of a data analysis method is to use ROC {Receiver 
Operating Characteristics). Such diagrams present the detection probability of a given signal (scaled at a particular 
SNR) versus the false alarm rate. 

Section II summarizes the general framework of this study, specifying the network detection model used in all the 
simulations presented in this paper, and the tools used. The single detector performances - detection efficiency and 
timing resolution - are briefly recalled in Section III as they are the basis for the following investigations. Then, 
the network coincidence analysis is studied in Section IV. ROC studies for networks from two to six detectors 
comprising the large interferometers currently being developed, GEO600 [1], the two 4-km LIGO detectors and the 
2-km interferometer located in Hanford [2], TAMA300 [3], Virgo [4] and finally the foreseen ACIGA project [5] - are 
presented. To decide whether or not a coincidence is valid, two different compatibility conditions are considered. The 
first one, called 'loose' condition, does not require the knowledge of the source position in the sky while the second 
one, the 'tight' condition, does. 

Section V deals with the coherent analysis. First, the derivation of the likelihood statistics follows closely the 
analysis of Ref. [20]; then, the corresponding ROC for coherent search of GW bursts are presented. As coherent 
filtering requires the knowledge of the source location in the sky, a set of coherent filters must be used in parallel to 
cover the full celestial sphere. Therefore, Section VI estimates the number of such templates needed for a complete 
tiling of the sky, given a prescription on the maximal loss of SNR allowed. Knowing allows one to tabulate the false 
alarm rate of the coherent ROC. Finally, Section VII compares the two network data analysis methods considered in 
this article. 



II. HYPOTHESIS AND NOTATIONS 
A. Interferometer response to a GW 

Measuring the strength of a signal with respect to the background noise is not the only information needed to 

estimate how well a GW may be detected in an interferometric detector. As the angular pattern of an antenna is 
not uniform, it is also necessary to take into account the location of the source in the sky. The result h{t) of the 
interaction between the wave and the instrument is a linear combination of the two GW polarizations h+ and hx [19]: 
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h{t) = F+{t)h+{t) + F^{t)h^{t) 



(2.1) 



The two weighting factors F+ and Fx are called beam pattern functions whose values are between -1 and 1. They 
depend on many parameters which may be roughly classified into three sets. 

(51) The detector coordinates (longitude, latitude, orientation with respect to the local North-South direction). 

(52) The source location in the sky, given for instance by the celestial sphere coordinates (the right ascension a 
and the declination S) and the local sidereal time which takes into account the Earth proper rotational motion. 

(53) A vector of physical parameters describing the time evolution of the GW signal. Some may be estimated 
at the output of the data analysis procedure, but this always requires adding some hypothesis on the signal. The 
polarization angle tp is also included in this set. 

In fact, the signal dependence on tp can be explicitly extracted. For an interferometer, one has: 

\F^{t)) \^-sin2^ cos2V'J \b{t)) 

where a{t) and b{t) only depend on the sets (SI) and (S2) [13,25], and with x being the angle between the two arms 
of the interferometer. 



B. Interferometer network and source modeling 

Monte-Carlo simulations [14] are used to compare the detection performances of coincident and coherent data anal- 
ysis methods in various network configurations. As this paper aims at studying the consequences of the interferometer 
location on Earth rather than the effect of the current or foreseen differences in their sensitivities to GW, the network 
model uses the simplifying assumption of identical detection performances. In this way, all interferometers contribute 
equally to the network. 

Yet, in Section IV C, a difference in the interferometer sensitivities is introduced in a complementary study, inves- 
tigating the case of non-identical detectors. The 'LIGO network' (made of three interferometers, the two 4-km in 
Hanford and Livingston, and the 2-km in Hanford) is well-suited to this work: as the GW sensitivity should scale 
with the arm length, the 2-km detector should ultimately be half as sensitive as the two other LIGO interferometers. 
As this computation shows an important loss of efficiency induced by this difference in sensitivity, the Hanford 2-km 
detector is not considered elsewhere in this article. 

The (SI) parameters of the detectors are chosen to match the already existing or planned instruments; as the local 
orientation of ACIGA is not yet defined, it has been optimized to maximize the detection efficiency in the full network 
of interferometers [14,24]. The P interferometers, labeled in the following by the index i, are assumed to have many 
features in common: the interaction with a GW signal - as defined in Section II A above -, the sampling frequency 
/samp, and the noise characteristics. All noises arc taken to be Gaussian, white and uncorrelated with the same RMS: 
(7* = cr. Finally, the interferometers are assumed to be properly synchronized. 

Any correlation between a filtering function s{t) and the output a;'(t) of the i-th interferometer at time t (the 
sampling time at the origin of the analyzing windows) is represented in the time domain by the following quantity: 

{s\x'){t) = J2s(j^) xx^(t+-^) (2.3) 

j^_Q \ /samp / \ /samp / 

with N being the filtering window size. 

The GW signal s{t) is assumed to be a Gaussian peak of 'half-width' iv = 1 ms: s{t) oc exp(— t^/2a;^). Such 
pulse-like shapes are characteristic of the most common GW burst waveforms simulated numerically: see e.g. [26,27] 
for the case of supernova signals. Its amplitude with respect to the background noise is monitored by its optimal 
SNR pmax [14]: the average value of the filter output, computed in a noisy background with both the Wiener filtering 
method and a detector optimally oriented. In the following, wc mainly focus on 'weak' signals for which detection 
problems are likely. We also assume that the sources are uniformly distributed over the sky, with a random timing. 
Finally, matched filtering is used to simulate the detection process. 
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C. Receiving Operator Characteristics (ROC) 



A convenient way to estimate the capability of a given filter to detect a particular GW signal - whose amplitude 
is fixed according to a chosen value of the maximum SNR pmax is the ROC, which presents on a single plot the 
detection efficiency e versus the filter false alarm rate r obtained by varying the threshold r]. Following the prescription 
of Ref. [18,28], we consider in the rest of this paper a false alarm rate per bin: 



"^""^ ~ N X iVMC ^ ' 

with A^FA being the total number of false alarms, A^mc the number of Monte-Carlo simulations, and N the size of the 
vector of data analyzed at each simulation loop. 



D. Coincidence analysis 



For the coincidence analysis framework presented in this paper, the event compatibility is tested by comparing the 
delays between the triggered events for every pair of detectors. In case of a real GW signal (detected in two detectors 
and at times and V respectively), the time difference is related to the source position in the sky. Let n be 
the unit vector pointing from the Earth center to the GW source location. One has: 



Af\n) = -f = (2.5) 

c 

if the filters have triggered on the GW signal in both detectors. Neglecting timing errors, such equation defines a 
circle in the sky on which the source is located. 

Two compatibility tests can be set from Eq. (2.5) for coincidence analysis: the first one - the 'loose' test - does not 
assume that the source location in the sky is known while the other - the 'tight' test uses this additional information. 
The former case allows one to survey the whole sky with one single analysis, but at the price of a lower efficiency, 
whereas the latter can reject more false coincidence events with a more stringent compatibility condition. 



1. Loose compatibility 



Let us first consider a full sky search without any knowledge on the GW source location in the sky. The timing 
delay At'-' between detectors and must obey the following inequality: 



< N N ^ ^ (2.6) 

c 

with the term on the right side of the inequality being the light time travel between the two detectors. Table I 
shows the maximum delays between all pairs of interferometers. The largest distance is between LIGO Livingston 
and ACIGA, about 42 ms. 



LIGO Hanford LIGO Livingston GEO600 TAMA300 ACIGA 


Virgo 27.20 26.39 3.20 


29.56 


37.06 


LIGO Hanford 10.00 25.01 


24.86 


39.33 


LIGO Livingston 25.04 


32.24 


41.68 


GEO600 


27.80 


37.46 


TAMA300 




24.58 



Table I: Maximum time delays (in ms) At^^^^ between pairs of interferometers. 
To take into account the statistical uncertainty on the timing locations, an error must be associated with the delay 
At^^ . For each interferometer, a single detection error is computed by using the relation between At^MS and the 
maximum filter output p shown on Figure 4 - see Section IIIB for more details. The two errors, assumed to be 
independent, are then quadratically summed. 
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The validity of the coincidence is tested by requiring \Afy\ to be smaUer than At^^^^ + ?]Jf°i5°g x AtJ^j^jg with 
^t?ming being a tunable positive parameter. Finally, multi-fold coincidences require that all participating pairs of 
interferometers are valid. Despite the apparent weakness of condition (2.6), taking into account all the delays available 
between detectors (a redundant set of informations) should nevertheless strongly cut false alarm events in the network 
considered. 

2. Tight compatibility 

If the source location in the sky is now a priori known, the compatibility condition can be tightened as the true 
delays At^-J^g between any pair interferometers are directly computed from Eq. (2.5). In this favorable case, the test 
requires the residual delay |At'^ — At^^^^j to be smaller than r/^l^ng ^ ^*rmS' with Jj^J^^g to be tuned as well to 
maximize the detection efficiency at a given false alarm rate. 

E. Simulation procedures 

Consecutive outputs of any burst search filter are highly correlated as the input data segments strongly overlap. 

Thus, algorithm outputs cannot be considered as statistically independent realizations of the same random variable: 
one often finds chisters of consecutive data exceeding a given threshold [10] which, in the case of a real GW burst, all 
correspond to the same signal. So, one has to redefine the event concept by counting only one single trigger when a 
consecutive set of ffiter output values are above the threshold. The two next paragraphs aim at giving some details 
on the Monte-Carlo simulation procedure for both coincidence and coherent analysis. Indeed, the latter method is 
more straightforward as all chunks of data are 'merged' in a precised way. 

1. Coincidences 

In the coincidence analysis simulations, an event is defined as a triplet of data: the maximum filter output, its 
associated time and the label of the interferometer in which it occurred. The coincidence ROC shown in this article 
have been constructed by using two different simulation steps: one computing the false alarm rate Tnorm, the other 
estimating the detection efficiency in the various network configurations. In both cases, the compatibility between 
alarms in different detectors is tested according to the prescriptions given in Section II D above. Results for loose and 
tight coincidence analyzes are presented in Section IV. 

Concerning false alarms, a two-step process is used in order to limit the computing time needed for simulation. First, 
the rate of (clustered) false alarms as a function of the triggering threshold is computed for the single interferometer 
case. Figure 1 shows in this case the evolution of Tnorm versus the threshold in a given window and some horizontal 
dashed lines translate Tnorm into more convenient values. Given the value of the threshold rj, this curve is used to 
generate random false alarms in a particular detector with a uniform time distribution. Finally, coincidences are 
searched in the lists of events associated with the different detectors in the network. 

Another point worth being mentioned is that the main effect of the alarm clustering procedure is to strongly reduce 
the false alarm rate Tnorm with respect to its estimator based on the assumption that consecutive filter outputs are 
independent. As shown in Figure 2, the ratio between the latter quantity and Tnorm is always above 10 and is indeed 
equal to the mean size of the false alarm clusters. 

To compute the detection efficiency, the first part of the simulation changes. As the maximum delay between 
two existing interferometers is 41.7 ms and as millisecond bursts are considered, synchronized data windows of N = 
1024 bins (corresponding to 51.2 ms for the sampling frequency /samp = 20 kHz) are enough to contain all the signal 
components after interaction with the detectors. Then, events are searched for in the data chunks from different 
detectors and the tests of coincidence compatibility between alarms are performed as in the false alarm case. 

As false alarm rates are kept low, the noise realizations are in this case required not to produce any false alarm. 
Of course, this bias is important only for large values of Tnorm, let say above 10~^. Below, the probability to have a 
false alarm in N data is under 1%, and so the noise bias does not play any significant role. In this way, one is sure 
that all triggers are due to some signal components and that coincidence efficiencies are not affected by false alarm 
contributions. 
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2. Coherent analysis 



The coherent analysis simulation is much simpler as the interferometer filter outputs are merged in a single data 
flow. Assuming a fixed analysis window size (N = 1024 for instance), false alarm rates and detection efficiencies are 
simply equal to the ratio of the number of simulations with outputs exceeding the threshold to the total number of 
simulations. The only assumption made is that errors in the source location are negligible, so that the relative shifts 
applied to synchronize the interferometer data segments are exact. 

III. SINGLE INTERFEROMETER STUDY 

This section summarizes the performances of a single interferometer, both in term of detection efficiency and of 
timing resolution, providing the necessary inputs for the coincidence studies. 

A. ROC 

Figure 3 collects some ROC for a single interferometer - the detection efficiency e versus the normalized false alarm 
rate Tnorm- Five different values of the optimal SNR Pmax are considered, ranging from 5 to 15. Like for all the similar 
plots of the article, some particular values of the normalized false alarm rate converted in more practical units are 
represented by vertical lines: from 1/year to 1/hour. This false alarm range should cover all interferometer operating 
configurations. 

Assuming pmax = 10 and an interferometer optimally orientated would lead to a detection efficiency very close to 

100% in the whole false alarm range. Unfortunately, because of the non-uniform antenna pattern, the detected signal 
amplitude is strongly reduced on average, and so the probability of detection. Thus, for intermediate values of Pmaxj 
the detection is not likely in a single detector, provided that the trigger threshold is kept high enough to have a small 
false alarm rate. Let us consider for instance the curve corresponding to pmax = 10, a typical value one can expect 
for a supernova at the Galactic center [26,27]. In order to reach a 50% efficiency, the detector must be run at a false 
alarm rate of roughly 1/second! 

B. Timing performances 

As all compatibility tests for network data analysis methods are based on time delays between the different interfer- 
ometer candidates, the timing resolution of a filter is another important quantity. Like for the detection problem, the 
matched filter appears to have the best resolution, as shown in [28] where optimal and sub-optimal filtering methods 
are compared. A first study of the Wiener filtering for a Gaussian signal in [14] showed that the timing resolution 
(i.e. the RMS of the difference At between the timing of the maximum filter output and the real GW timing) could 
be simply parameterized: 

AtHMS« 0.15ms (3.1) 

with an excellent agreement between the fit and the real RMS as soon as p > 6. 

Here, we extend this work by studying the evolution of A%ms on a larger range of p: from (no signal) to a very 
large value - see Figure 4. This result will be used later to validate coincidences in Section IV. 

As we fociis in this article on a Gaussian peak with width u = 1 ms lasting around 6 ms in total, a window of 
N = 512 has been chosen to compute the evolution of AfRMS versus p. Indeed, it corresponds to 25.6 ms for a 20 kHz 
sampling frequency, a duration large enough to include the whole signal. Choosing a much larger value for N is not 
suitable as for a negligible GW signal completely dominated by the noise, At is uniformly distributed in the analysis 
window. The timing error RMS would then grow 'artificially' with N. Consequently, the compatibility condition 
would be more easily satisfied, leading thus to an increase of the false alarm rate. On the other hand, for very large 
values of p, one expect AtpMS to scale like 1/^^, i.e. a slower variation than Eq. (3.1). Ultimately, the timing 
resolution would be limited by the sampling frequency. 
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IV. COINCIDENCES 



In this section, different configurations of interferometric detector networks are studied: the 3-interferonieter network 

Virgo-LIGO, the LIGO 3-intcrfcromctcr network inchiding the Hanford 2- km detector, and the full network made of 
the six first generation interferometers. They allow one to predict the performances that could be achieved in detecting 
GW bursts in coincidence in the future. This network detection model could easily be updated when final relative 
sensitivities are known, or when new second generation detectors appear. In these studies, the loose compatibility 
test is always used, apart in the last paragraph IV D where a summary of tight coincidence performances is presented. 

A. Virgo-LIGO network 

First, we consider the three-interferometer network Virgo-LIGO. Studying its efficiency is important for two main 
reasons: it includes the detectors with the best foreseen sensitivities and a threefold detection is the minimum number 
of coincidences required to estimate the source location in the sky. 

1. Two-interferometer coincidences 

The simplest network is made of two interferometers. Its performances depend a lot of the particular configuration 
considered, as shown in the following. For the timing compatibility condition, three different values of tiI°^^^ have 
been tested in simulations: 1, 2 and 3. It turns out that the best compromise between low false alarm rate and high 
detection efficiency is obtained with r/^^^^g = 1, value used in the following for all loose coincidence tests. 

Figure 5 compares the ROC computed for the three pairs of detectors with pmax = 10. The first point to notice 
is that the efficiency never reaches 60%, even at very high false alarm rates: two interferometers are not enough to 
guarantee a likely detection of such bursts. Then, the configuration associating the two LIGO detectors shows clearly 
better performances than the two others made of Virgo plus one LIGO. Two reasons explain these differences: 

• The two LIGO interferometers have been built in order to maximize the correlation between their antenna 
patterns, increasing the coincidence efficiency. This is the dominant effect - see the next paragraph. 

• The LIGO detectors are close with respect to Virgo - see Table I in Section II D - and less random coincidences 
are allowed in the compatibility window. So, the false alarm rate is shifted to the left thanks to this effect. 

This feature remains true if one compares all pairs of interferometers chosen among the full network of six detectors: 
the two LIGO configuration ROC is clearly better than any other. Even the Virgo-GEO600 pair - the two closest 
instruments - cannot compete: for a given threshold, the false alarm rate is lower, but also the efficiency as the 
angular patterns do not overlap well. 

Figure 6 shows how the two LIGO detection efficiency evolves for different values of the optimal SNR, between 5 
and 15. The larger Pmax, the better the efficiency, but the improvement remains limited: even for Pmax = 15, the 
detection probability is only around 50% for Tnorm = 1/ hour. So, a two detector network appears to be not sufficient. 
Yet, a last point to be mentioned is that the efficiency decreases more slowly with the false alarm rate than for the 
single detector case - see Figure 3 for comparison. The next section will clearly show the interest of this behaviour. 

2. Coincidence strategy comparison 

Figure 7 compares for pmax = 10 all the possible coincidence strategies in the Virgo-LIGO network: single detector, 
coincidences in the two LIGO interferometers (the best pair of detectors) , twofold coincidences (at least two detections 
among three) and finally full coincidences. The twofold coincidence strategy is clearly the best: its ROC is above the 
other ones in the full range of false alarms covered by the graph. Yet, it does not show very high efficiencies: only 
a bit more than 40% for 1 false alarm per hour, and around 25% at the level of 1 per year. In addition, threefold 
coincidences are very inefficient, and these results are indeed similar for all triplets of interferometers. The networks 
with the best ROC all include the two LIGO detectors, but their efficiencies depend weakly on the location of the third 
one on Earth: replacing Virgo by GEO600 gives a slightly better result, while using TAMA300 or ACIGA decreases 
a bit the efficiency at given false alarm rates. Therefore, a three interferometer network does not appear large enough 
to reach high detection efficiencies. 
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Two other interesting points can be extracted from the plots. First, as soon as the false alarm rate becomes small 
enough to be realistic (say around Tnorm = 1/ few seconds), the ROC corresponding to coincidences between the two 
LIGO detectors crosses the single interferometer ROC and shows larger efficiencies at fixed false alarm rate. This 
crossing is due to the fact that the large decrease in the false alarm rate due to the compatibility condition required 
for a LIGO coincidence is stronger than the corresponding loss in efficiency caused by the twofold detection. This 
demonstrates that searching GW bursts in one detector is not only discarded by the small confidence level one can 
associate to such events, but also because this method is simply less efficient than others. 

Second, for smaller false alarm rates, the ROC for twofold coincidences and for the two LIGO get closer. This is due 
to the fact that the two LIGO antenna patterns are close one each other and quite different from the Virgo pattern. 
Therefore, when the threshold increases, coincidences between Virgo and one LIGO are more strongly suppressed 
than for the LIGO pair. Yet, the two curves are widely separated in the whole range of false alarm: adding Virgo 
allows one to improve significantly the efficiency with respect to the two LIGO detectors alone more than 30% on 
a relative scale. So, going from two to three interferometers in the network is a clear improvement. 

B. Full network of six interferometers 

As the previous section pointed out that a three interferometer network is not promising enough, we compute in 
this section the detection efficiency for the full network of six detectors. Operating such configuration in a near future 
should be a clear goal of the worldwide GW community. Indeed, this network reaches quite promising detection 
efficiencies, as Figure 8 shows. The ROC, computed again for pmax = 10, correspond to coincidence strategies in 
which a minimal number of detections is required: from two (top curve) to six (bottom curve). In this configuration, 
twofold coincidences arc quite likely: more than 80% for 1 false alarm per hour and still about 60% for 1/year. In 
addition, threefold coincidences appear possible: in this case, the efficiency is around 60% for 1 false alarm per day. 
On the other hand, higher multiplicity coincidences are less and less efficient. 

Table II summarizes the loose coincidence results previously presented. It collects detection efficiencies from various 
strategies, sampled at representative false alarm rates. 



Configuration 


Tnorm = 1/year 


Tnorm = 1/week 


Tnorm = 1/day 


Tnorm = 1/hoUr 


Single detector 


16% 


20% 


23% 


30% 


At least 2/3 (Virgo-LIGO network) 


27% 


34% 


37% 


44% 


3/3 (Virgo-LIGO network) 


10% 


13% 


15% 


20% 


At least 2/6 


57% 


69% 


74% 


83% 


At least 3/6 


47% 


52% 


57% 


65% 


At least 4/6 


23% 


31% 


35% 


42% 



Table II: Loose coincidence efficiency comparison for pmax = 10 



C. The LIGO network 

The LIGO system actually consists of three interferometers: the two 4-km detectors in Hanford and Livingston, 
and the Hanford 2-km interferometer located in the same vacuum tube as the larger instrument. Assuming the most 
optimistic situation in which the dominant noises of these two neighbor detectors are independent, their close locations 
significantly reduce the number of random false alarms between them. On the other hand, the difference in the arm 
lengths reduces the 2-km detector sensitivity by a factor two. Therefore, it is interesting to see how these two effects 
balance and ROC are well-suited for such a study. 

Figure 9 compares the full LIGO network with the Virgo-LIGO (4-km) network. The top plot presents ROC for 
pairs of interferometers: the two LIGO 4-km detectors, Virgo and each of the LIGO 4-km interferometers, and finally 
each of the LIGO 4-km interferometers with the Hanford 2-km detector. The bottom plot compares the coincidence 
strategies involving the three detectors in each network: twofold coincidence (at least two detections among three) 
and the full coincidence. 

These graphs show that the performance of the full LIGO network is worse than the Virgo-LIGO (4 km) network: 
detections efficiencies at given false alarm rate are better in the latter case. The reduction factor of the LIGO Hanford 
2-km detector sensitivity plays a more important role than the gain in false alarm rate provided by the coincident 
locations of the Hanford interferometers. Conversely, this is a strong indication that adding in a network detectors 
less sensitive than others will only give limited improvements in detection efficiency. 
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D. Tight coincidences 



To conclude this section dealing with coincidence detection, Figure 10 compares ROC computed for different tight 

coincidence strategies: coincidences between the two 4-km LIGO detectors, twofold and threefold coincidences in both 
the Virgo-LIGO network and the full network of six interferometers. As the source location is assumed to be known 
here, such curves can be directly compared with the coherent analysis results presented in the next section. 

As for the loose coincidence case, a tuning of the parameter r?*i'^*ig has been performed; values ranging between 1 

and 5 have been used, with the choice TyJ-^^g = 3 giving the best ROC. Comparing Figure 10 with Figures 5 to 8 
- corresponding to the loose compatibility criterion - shows that the tight coincidence ROC present larger relative 
efficiencies. Yet, the improvement remains limited, from 10 to 20% in relative. Therefore, the main limitation of the 
coincidence ROC appears to be the coincidence algorithm itself, handling only binary informations (signal present or 
absent) in the different detectors with a fixed threshold. 

V. COHERENT ANALYSIS 

Compared to coincidence searches, coherent analysis methods use a more complete set of informations coming from 
the different components of the network; thus, their detection efficiency is larger. On the other hand, using them in 
real analysis requires an additional hypothesis on the source location, which allows one to properly shift the various 
detector outputs to synchronize them. In the most general case, the source position in the sky, monitored by two 
angles, is unknown and thus must be added to the set of unknown parameters describing the GW signal. So, many 
templates must be run in parallel to ensure an efficient coverage of the sky, each of them focusing on a particular 
area. Consequently, any comparison between coincidence and coherent methods must take into account this fact, 
leading in particular to a renormalization of the false alarm rate in the case of loose coincidences see Section VII 
for more details. Yet, if the source location is already known (e.g. from informations given by detectors sensitive to 
other radiations), this restriction is lifted and one coherent algorithm is enough. In this case, its performances can be 
directly compared with the tight coincidence scenario. 

To cope with this requirement, the study of coherent analysis methods is performed in two steps. In this section, 
ROC are computed for a single coherent algorithm, assuming a perfect knowledge of the source location in the sky. 
Then, the number of filters needed to cover the sky is estimated in Section VI, using the formalism of Ref. [23]. 

A. Derivation of the statistics from the likelihood ratio 

In the following, the derivation of the coherent statistics based on the likelihood ratio is briefly recalled. Assuming 
known the source sky location, the template to be used in the i-th detector takes the form 

s\t) = Fi X X s°{t - Sf) (5.1) 

beam pattern term « 1/distance 

where s° is the generic template shifted by the time delay dt^ and scaled by Fi, the factor giving the quality of the 
interaction between the antenna and the GW. Without loss of generality, one can assume that (s" | s°)(0) = 1. 

Let us consider first the search of a known signal in a single detector output. The most efficient method is in 
this case the Wiener filter. Its expression naturally arises in the framework of the likelihood ratio, defined as the 
conditional probability to have a particular set of data assuming that the signal is present - see e.g. the corresponding 
discussion in [18]. 

By using the same method, one can define a global likelihood ratio for a set of detectors. With the hypothesis made 
on the interferometer noises (Gaussianity and independence), the logarithm of the 'network' likelihood ratio InA is 
computed by simply summing the corresponding contributions for the P single detectors D*: 



InA 

Its final expression can be found in Ref. [18,20]: 



^InA^ 
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i=l 



i=l 



(5.2) 



which is a second order polynomial function in the unknown K. Still following [20], maximizing over this variable 
and extracting the polarization angle V' from the beam pattern functions gives the following expression: 



In A 



(cos2V'A + sin2V'B) .' 
II cos2V'^ + sin2V'-B|| 



(5.3) 



with: 



a- 



and f (to) 



for 



1,...,P 



The vector 7 can be expanded in the following way: 



7 



TaA + TbB + j± with j± ± A and j± ± B 



(5.4) 



Therefore, only the two first terms of this sum contribute to the likelihood ratio. One can now compute the value 
of tp which maximizes it. By applying the Schwarz inequality, one gets the following upper bound 



(In A) 



< 



TaA 



TbB 



(5.5) 



which is reached by choosing 2^' such as 



cos(2'^) 



and sin(2^) 



The maximum of the likelihood ratio statistics is thus proportional to the square of the norm of the orthogonal 
projection of 7 on the subspace generated by the couple (^A, B^ . Computing the values of Ta and Tb allows one to 
give a compact expression of this new statistics, denoted A in the following: 



A 



B] A - 



A] B 



A^B'^ - [a. By 



(5.6) 



with A ' 



and B 



B 



This last expression shows that all the former calculations are meaningless if A and B are parallel ~ indeed, 
invalidating Eq. (5.4). Fortunately, this critical situation is unlikely in the sky, as shown in Figure 11 where the sky 
map of I cos^^bI = \A.B\ / (AB) is represented as a function of the celestial sphere coordinates for the Virgo-LIGO 
network and for the full set of six detectors. 

For the first network, the values cos 9 ab = ±1 can be reached only in a very small area of the sky while in the 
second case, the absolute value of the cosine remains below 0.6 in any direction of the sky. This difference is simply 
due to the fact that the vectors A and B have grown from 3 to 6 components and are therefore less likely to be 
collinear. 
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B. Statistical behavior of A 



To use the estimator A for data analysis purpose, one can first study its statistical properties under the hypothesis 
of noise only. To do this, one has to rewrite Eq. (5.6) on a different way. First, introducing an orthogonal basis in 
the plane containing A and B by defining the two following vectors [20] : 



A can then be rewritten: 



A = 



^ A B ^ ^ A B 
u=- + - and v = - - - 



^2b2||^||2 ||^||2 



4 A^B"^ - [A.B 



^ \ 2 

^'W\\. 



(5.7) 



Coefficient only depending on tlie source slcy location 

The second term of the previous equation is the reduced statistics used in the following: 



A, 



reduced — 



7- 



From the definitions of the vectors 7, u and v, it clearly follows that the two Gaussian variables T-j]^ and T-jifij 



• 11^711 



uncorrelated. Indeed, the distribution of Areduced is close to a variable with two degrees of freedom, independently 
of the source location in the sky and of the particular network considered. 



C. Coherent data analysis ROC 

Figures 12 and 13 present two examples of coherent ROC, for the three detector network Virgo-LIGO and for the 

full set of six interferometers. For the two networks considered here, the coherent analysis ROC are clearly above 
all coincidence ROC. This is mostly due to the more complete management of data in the coherent method case. 
Moreover, comparing with the loose coincidence case, the coherent approach benefits in addition from the fact that 
the source location is known. Section VII summarizes the comparison of both network data analysis approaches. 

One can also note that going from three to six interferometers strongly increases the detection probabilities; the 
differences are more significant than for the coincidence case. Indeed, for Pmax = 10; the efficiency remains higher 
than 97% in the whole range of false alarm rates per bin; for Pmax = 7.5, the detection efficiency is at least 80%. For 
T = 1/hour and Pmax — 5, one has still e = 50%. 



D. Coherent data analysis timing accuracy 

As for the single detector case, the timing accuracy of the coherent method can be easily studied - estimating the 
timing accuracy for coincidences is not as straightforward. Figure 14 shows the evolution of the timing error At^MS 

1 /2 

as a function of (Areduced) for the two examples of networks considered in this section: Virgo-LIGO and the full set 
of six interferometers. As for the single detector case, the precision goes well below the signal half-width to - taken 
equal to 1 ms here. Taking the square-root of the (quadratic) coherent statistics is mandatory in order to have a 
quantity scaling with the optimal SNR Pmax- Due to the 'universality' of Areduced, the two curves presented on the 
plot overlap perfectly. 

1/2 

One can also try to connect (Areduced) and Pmax, at least in average. For pmax > 3, linear fits give: 



1 /2 

(Areduced) ~ 0.67 X Prnax + 0.89 for the Virgo-LIGO network 

1 /2 

(Areduced) w 1.06 X Pmax + 0.35 for the six- interferometer network 
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Inverting those cqnations allows one to roughly link (Aj-cduccd)^^^ to a value of the optimal SNR pmax- Of course, the 

1 /2 

larger the network, the higher the mean value of (Areduced) at fixed Pmax- 

The proper way to estimate the timing accuracy improvement provided by a network coherent analysis with respect 
to the single interferometer case is to compare the timing error RMS for a given optimal SNR Pmax- Table HI below 
shows the timing performances of the different configurations for three values of Pmax- 5, 7.5 and 10 respectively. 
Only events exceeding the threshold tuned at a false alarm rate of 1 /hour are included in the computation. 





Single detector 


Virgo-LIGO network 

(coherent analysis) 


Full network 
(coherent analysis) 


5 


0.44 ms 


0.29 ms 


0.25 ms 


7.5 


0.29 ms 


0.23 ms 


0.18 ms 


10 


0.24 ms 


0.19 ms 


0.14 ms 



Table III: Timing performance comparison for different values of Pmax- 
As expected, the coherent analysis improves also the timing precision, especially at low optimal SNR. GW events 
do not only trigger more often; they are also more precisely located. Finally, Table IV and V show how the coherent 
timing performances at given pmax evolve when the false alarm rate is reduced. As the thresholds increase, the quality 
of the selected sample improves; yet, the precision in locating the GW signal peaks does not change significantly. 



False alarm rate 


Pmax — 5 


Pmax — 7.5 


Pmax = 10 


1 /hour 


0.29 ms 


0.23 ms 


0.19 ms 


1/day 


0.25 ms 


0.21 ms 


0.18 ms 


1 /week 


0.24 ms 


0.20 ms 


0.17 ms 



Table IV: Coherent timing accuracy for decreasing false alarm rates in the Virgo-LIGO network. 



False alarm rate 


Pmax — 5 


Pmax — 7.5 


Pmax = 10 


1 /hour 


0.25 ms 


0.18 ms 


0.14 ms 


1/day 


0.24 ms 


0.18 ms 


0.14 ms 


1 / week 


0.23 ms 


0.18 ms 


0.14 ms 



Table V: Coherent timing accuracy for decreasing false alarm rates in the full network. 

VI. MATCHED FILTERING OF THE CELESTIAL SPHERE FOR A COHERENT ANALYSIS WITH A 

NETWORK OF INTERFEROMETERS 

The last step of the coherent search - when the source location is a priori unknown - consists in estimating the 
number of filters 9^ needed for the sky coverage. To do this, the most efficient way is to use the method first defined 
in Ref. [23] for the in-spiral binary case, and then extended for coherent analysis of Newtonian chirp binary signals 
in a network up to three interferometers [20]. Here, we still use a Gaussian peak of width uj as 'generic' GW burst 
signals. 

The main difference with the matched filtering case is that there is no more symmetry between the interferometer 
data and the template. Quantifying the separation between two close filters is thus not easy. For instance, let us 
consider the case of the polarization angle i]r. as shown in section V A. there exists an analytical way to maximize A 
over ^/'template while ^/'signal remains 'hidden' in the noisy data. The solution proposed by Ref. [20] is to choose some 
values for Vsignai - and also for the binary orbit inclination in that case - and to estimate the number of templates 
^Jl for these different configurations. Numerically, it is found that 9T does not change by more than a factor 3 in the 
range of parameters tested. In this paper, a different path is followed: the logarithm of the likelihood ratio is first 
averaged over xp, which allows one to focus only on the sky angular dependence of the beam pattern function. 

Before presenting the calculation one can remark that, as first pointed out in [20], the loss in SNR caused by a 
mismatch in the source direction is mainly due to the corresponding wrong time delays. 

A. Ambiguity function and metric in the celestial coordinates 

Assuming a mismatch St' between the Gaussian peak template and the GW signal - both of characteristic width 
a; -, a straightforward calculation of the correlation gives: 
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{s'\x') = {F'f exp 



6V_ 
2w 



(6.1) 



The exponential term reduces the maximal correlation and the signal width lo provides a timescale with which the 
error 5t^ is compared. Of course, Eq. (6.1) would be meaningless in the case of a single detector, as a simple 
opposite time-shift of the template would allow one to recover the full SNR. In coherent analysis, time shifts cannot 
be optimized separately for each interferometer: therefore, a wrong match of the detection leads to unavoidable losses 
in SNR. The logarithm of the likelihood ratio averaged on the polarization angle is thus equal to: 



InA = K^Y. i^'y 



exp 



V2w 



with = — r 



(6.2) 



To 'transform' In A into an ambiguity function giving the relative mean loss in SNR due to the direction mismatch, 
one chooses K such as lnA|(5ti=o) = 1) which is achieved with: 



K = 



Ef=r {^'f 



A Taylor expansion A around ((5t* = 0)i=i....^p at the second order gives the quadratic approximation of the 
ambiguity function, assumed to be valid provided that the allowed losses of SNR remain small. 



1 Etr {^'f m 



with 



p 

E 

i=l 



(6.3) 



Then, one has to replace the by their expressions in term of the two angular variables locating the source 
direction in the sky: the right ascension a and the sine of the declination X = sm5. Let Q. represent the center of 
Earth and n be the unit vector radiating from it in the source direction. One has 



5f = — 5n. no' 
c 



with c being the speed of light and 6n the error in the direction of the source location - note that from Eq. (6.2) the 
sign convention of 6t^ does not matter. The computation of 6n is straightforward: 



X 



Sn = 



-Vl - sina --^=^ 
VI - cos a 







1 



da 
dX 



The 3x2 matrix appearing in the previous equation will be designed as M in the following. In order to shorten the 
expressions appearing in the metric calculation, one can also introduce a 3 x 3 matrix F defined as follows: 



with 



The r matrix contains all the network characteristics. Equation (6.3) can thus be rewritten: 



1 



Sn] .T.Sn 



1 



4 c^ ip^ 



{da dX) . (*Mr: 



da 
dX 



(6.4) 
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The tiling metric g is thus: 



_ G 
^ ~ 4a;2 c2 (^2 

The exact expressions for the coefficients of the (symmetrical) 2x2 matrix G are given in Appendix A, so as its 
determinant Ac- 



B. Number of templates needed for the various network configurations 

Calculating the metric allows one to estimate the mimbcr of templates '^20 needed to cover the whole sky for a 
particular network of interferometers, given the maximal allowed loss in SNR. This last quantity is usually written as 
1 - MM where MM is the 'Minimal Match' [23] - a conventional value is MM = 97%. 

Still following Ref. [23], is computed by integrating over the sky the square root of the metric determinant, 
multiplied by a scaling factor depending on the minimal match and on the parameter space {(a, X)}. One gets finally 



8a;2c2(l - MM) y[_.;.]x[-i;i] 



/ 

J f— 7r;7rl 



dadX 



(6.5) 



Prom this formula, one can note that the longer the signal, the smaller the number of templates. This last feature 
is due to the particular burst shape chosen in the paper and cannot be generalized to any GW signal (indeed pure 
sines behave in an exactly opposite way as the larger their number of cycles, the more they need to be accurately 
tracked in noisy data). Appendix B gives an estimation of OT for the different networks of existing interferometers. 

One can see that the smaller the network, the more the value of (Tl depends on the particular configuration. This is 
particularly true for the case of P = 2 detectors for which there is a factor higher than six between the extreme values. 
In this case, the number of templates does not only depend on the light-distance between the two interferometers but 
also on their respective orientations. For larger networks, the results are closer: the exact locations of the detectors 
appear less important, they look like more 'randomly' spaced on Earth. For the set of 6 interferometers, one has 
m ~ 5320. 



C. Extending the space peirameter 

One can also assume that the exact width of the GW signal is not known and thus that w is another parameter of 
the search. It is easy to check that the only change in Eq. (6.3) is the apparition of a new term reducing in addition 
the ambiguity function: 




where Suj is the error on the Gaussian peak width. As the width and the angular parameters are decoupled, estimating 
91 is straightforward: 

msn - . 'i^.^W2 ( I \^\ d-dx\ ( j %\ (6.6) 

64 c2 (1 - MM)3/2 \^ L ^2 J j yj^^^^^.^^^^^ J 

To measure the 'template cost' due to the addition of the third free parameter ui, one can for instance compute 
the ratio between the '3D' number of filters needed to fill both the w-range [wminiWmax] and the corresponding '2D' 
number for u) = Wmin fixed. As seen from Equations (6.5) and (6.6), this ratio does not depend on the network as the 
angular integrals simplify. 



^2D (Wmin) 



3\/3 



16\/1 - MM 



(6.7) 
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For MM = 97%, the numerical factor in front of the brackets is equal to 1.88 and this asymptotic value is quickly 
reached when the ratio i^min/t<-'max decreases. The number of templates only doubles when one goes from two to three 
parameters; thus, covering coherently the celestial sphere for a burst search over a wide range of durations is not 
too expensive. This number has to be compared with the overestimated value of computed by multiplying 
by the number of filters needed to cover the one-dimensional parameter space [wmin; Wmax]- As 91^; = 12 for the 
numerical data considered in this section [29], the saving in template number - and thus in CPU time - is at least a 
factor 6. 



VII. COMPARISON OF COINCIDENCES AND COHERENT DATA ANALYSIS METHODS 

As previously stated, comparing the ROC for coherent and coincident analysis leads to the clear conclusion that 
the former approach shows better performances. Indeed, methods based on a coherent use of the various datasets 
must give better results than coincidences, as the merging of informations provided by the different interferometers 
is more complete than a simple binary test (absence or presence of the signal in a given detector). 

A very strong assumption made for the study of the coherent analysis is that the source location is known, whereas 
no such hypothesis is necessary for loose coincidence detections. A priori, this additional information could be the 
main origin of the performance differences between the two network data analysis methods. Yet, the studies performed 
in this paper do not confirm this hypothesis of the dominant improvement factor. Indeed, tight coincidences have 
also been studied. In this case, the source location is assumed to be known, as for the coherent analysis, and so ROC 
are directly comparable. Yet, the performance gap between the two network algorithms remains wide. 

Requiring the knowledge of the source location for coherent filtering has also another consequence: to cover the full 
sky, many templates must be used in parallel. Therefore, the meaningful quantity is no more the false alarm rate per 
bin Tsingie of a given coherent filter, but rather the global false rate Tgiobah computed by taking into account the whole 
set of templates. As the parameter space grid is thin, false alarms between close filters are certainly correlated: if one 
filter triggers, some templates corresponding to neighbor locations should also exceed the threshold. Computing the 
correlation level is a complete work by itself; thus, in this article, we only estimated it roughly with a toy Monte-Carlo 
cheating on the precise location of the templates. Assuming 1 false alarm per hour and per template, the fraction of 
filters triggering simultaneously is Kcorrei ~ 7% for the Virgo-LIGO network and Kcorrei ~ 0.5% for the full set of 6 
interferometers. As 



with *n being the number of templates computed in the previous section. From the numerical results given in Appendix 
B, one can deduce that with a minimal match MM = 97% one has 



as the template number is OT ^ 5000 in both cases. 

Comparing Figures 3 to 8 on the one hand and Figures 12-13 on the other hand, clearly shows that even if the 
horizontal-axis of the coherent ROC are shifted by these values on the right, the corresponding detection probabilities 
remain clearly higher than for the coincidence methods. 

Indeed, Table VI summarizes the coincidence and (rescaled) coherent analysis efficiencies^ for different false alarm 
rates. Three scenarii are compared: 

• Tgiobai = Tsingic: no Correlation between templates; 

• tcorrci equal to the estimations presented above; 

• the worst (and unlikely) case, Kcorrei = 1^ maximum correlation. 



^The values tagged with an (*) in Table VI have been estimated by prolongating the ROC beyond the range of false alarm 
rates achieved by the numerical simulations. 



"^global '"^ '^single ^ '^correl ^ 



(7.1) 




350 X 

''single for the Virgo-LIGO network 

25 X 

"^single for the full network 



(7.2) 
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False alarm rate Tnorm 


1 / week 


1/day 


1 /hour 


Twofold loose coincidence 


34% 


37% 


44% 


Twofold tight coincidence 


37% 


41% 


47% 


Coherent analysis (no correlation) 


61% 


65% 


72% 


Coherent analysis (Kcorrei = 7%) 


46% (*) 


53% (*) 


60% 


Coherent analysis (Kcorrei = 1) 


43% (*) 


48% (*) 


54% (*) 



Virgo-LIGO network 



False alarm rate Tnorm 


1 / week 


1/day 


1 /hour 


Twofold loose coincidence 


69% 


74% 


83% 


Twofold tight coincidence 


76% 


82% 


87% 


Threefold loose coincidence 


52% 


57% 


65% 


Threefold tight coincidence 


61% 


64% 


72% 


Coherent analysis (no correlation) 


98% 


99% 


100% 


Coherent analysis (Kcorrei = 0.5%) 


98.1% (*) 


98.6% 


99.3% 


Coherent analysis (Kcorrei = 1) 


96.0% (*) 


96.7% (*) 


97.6% (*) 



Full network 



Table VI: Comparison between coincidence and (rescalcd) coherent analysis detection efficiencies at various false alarm rates. 

Finally, one can note that keeping Kcorrei constant when Tnorm decreases leads to an overestimation of the template 
correlations at smaller false alarm rates: the higher the threshold, the smaller the probability to have again the 
filter outputs triggering when the datasets are shifted one with respect to the other. Therefore, coherently analyzing 
data coming from the different detectors increases significantly in all cases the detection potential of interferometer 
networks. Moreover, the number of templates involved in such searches appears low enough to make these analyzes 
feasible with a small CPU farm. 

The main reason why coherent analysis appears so successful is certainly its capability to sum the signal contributions 
from the different interferometers regardless whether they individually trigger. A coherent detection can perfectly 
originate from outputs distributed in such a way that none of them triggers on coincidence strategies for thresholds 
adjusted to the same false alarm rate! On the other hand, coincidences always loose a significant fraction of the 
available information, which becomes more and more important as the network size increases: the larger the network, 
the more likely that a GW signal strong enough is above the background noise in some of these detectors. 

VIII. CONCLUSION 

Coincidence and coherent data analysis methods in networks of interferometric GW detectors are compared in this 
article thanks to a network model which allows one to compare these two approaches quantitatively, through ROC 
curves. Indeed, these graphs summarize well the behavior of a particular algorithm for a given GW signal over a wide 
range of false alarm rates. First, coincidence methods are studied in different networks from two to six interferometers. 
To select events, two different compatibility criteria arc defined. The first one, the loose test, does not require any 
assumption on the source location in the sky; therefore, using it allows one to search events in the whole celestial 
sphere, however with limited efficiency. 

From this study, it clearly appears that searching GW bursts in a single detector is not efficient at all. For what 
concerns two detector networks, the LIGO 4-km pair is the most efficient, due to their relative closeness and especially 
to their 'parallel' orientation. Yet, detection efficiencies remain limited for such reduced networks. Therefore, larger 
sets of detectors have to be considered. Adding Virgo to the two LIGO detectors shows a significant enhancement 
of the twofold detection efficiency. But higher-fold coincidences remain improbable, unless the network size increases 
significantly. So. the goal of a worldwide coincidence analysis should be to include as much interferometers as possible 
- indeed, six would not be too much! - in the network. 

A complementary study of the full LIGO network including the two 4-km interferometers and the 2-km detector in 
Hanford shows that differences in the sensitivity of network components have important consequences on the network 
performances: the factor two difference in the 2-km interferometer sensitivity is more important than the false alarm 
rate reduction due to the close location of the two Hanford detectors. Therefore, an efficient network should contain 
interferometers with sensitivities as close as possible. Conversely, adding less efficient instruments to a network is not 
productive. 

Finally, few ROC about tight coincidences are presented. As the source location is known, the compatibility test is 
more constraining, leading thus to an improvement of the ROC performances. Yet, the gain is small as coincidence 
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analyses arc limited by the loss of information due to the binary diagnostic made in each interferometer of the network. 
On the other hand, coherent analyses benefit from all detector outputs without a priori on the presence/absence of a 
GW signal in the data, and are so much more powerful. 

The difference in performances between coincidence and coherent analysis is significant in both networks considered 
here: Virgo-LIGO and the full network of six detectors. Coherent detection efficiencies remain very large even at 
small false alarm rate for Pmax = 10 - more than 95% at Tnorm = 1/week in the six detector network! On the other 
hand, weaker signals (say pmax = 5 or below) are not well seen: the final sensitivities of network components will be 
critical. Finally, the timing accuracy of coherent data analysis method is also considered; as expected, it is in average 
better than the single interferometer case, and improves with the network size. 

Another point worth being mentioned about coherent analysis methods is that they could be used even when the 
waveform is not accurately known, like for GW bursts in general. Indeed, the only assumption made in this paper is 
that the filtering algorithm is linear - cf. Eq. (2.3). For instance, one could also use some robust and efficient filters 
[6,28] which only depend on a single parameter, the analysis window size. 

The price to pay for the high performances of the coherent data analysis method is its complexity, especially with 
respect to the loose coincidence approach. But, at least in the case of burst signals, this does not appear to be a 
strong limitation: the number of templates needed to scan accurately the whole sky is quite small - at most a few 
thousands -, even including the signal width. In addition, the correlation between the templates is estimated to be 
below a few percent or even less. Therefore, the loss in performances induced by the increase of the global coherent 
false alarm rate with respect to the single template case is limited: coherent methods are better than coincidence 
searches. 

So, the main conclusion of this study is that one should not limit collaborative data analysis to the exchange of 
single interferometer events, especially for GW bursts. Otherwise, a large fraction of detection efficiency will be lost, 
which may be crucial for rare sources, like e.g. close supernovae. Using the full set of available data for GW signal 
search - in the largest possible network - should be an important goal of the worldwide GW data analysis community, 
at least in a mid-term perspective. 

APPENDIX A: COEFFICIENTS AND DETERMINANT OF THE METRIC MATRIX G 

With the notations defined in the core of the paper, the coefficients of the 2x2 symmetrical matrix G are: 
Gil = (1 - X^) [ Til sin^ a - Fia sin(2a) + F22 cos^ a ] 

+ Vl ~ ( r23Cosa — Fis sin a ) 



(Fii - F12) sin 2a 
^ — ^ — F12 cos 2a 



G\2 = G21 = X 

G22 = z ( Til cos^a + F22 sin^a + F12 sin 2a ) - (F13 cos a + F23 sin a) 

i- — X y/\ — X-^ 

To estimate the number of templates Vt, one needs to compute the determinant of G. Extensive calculations give: 

Ag = ( Fn F22 - Ffj ) X^ 

- 2X Vl - [(F22F13 - F12F23) cos a + (F11F23 - F12F13) sin a] 

+ (1 - X2) [ (F22F33 - Fy cos^a + (F11F33 - Ffa) sin^a + (F13F23 - F12 F33) sin 2a ] 

As shown by the previous formula, Ag is never singular, even in the directions X = ±.1. 



APPENDIX B: LIST OF THE NUMBERS OF TEMPLATES FOR THE DIFFERENT 

CONFIGURATIONS OF DETECTORS 

This appendix gives the estimated number of templates needed to cover the whole sky for each possible network of 
interferometers: from 2 detectors to the whole set of 6 antennas. These numbers are computed with w = 1 ms and 
MM = 97% and thus must be properly rescaled for different choices of these two coefficients - see Eq. (6.5). 

In the following tables compiling the results of the calculations, the interferometers are simply designed by a single 
letter: Virgo (V), LIGO Hanford (H) and LIGO Livingston (L), GEO600 (G), TAMA300 (T) and ACIGA (A). 
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FIG. 1. Normalized false alarm rate Tnorm as a function of the triggering threshold rj. This curve is computed for a window 

size N = 1024 bins, i.e. 51.2 ms at a 20 kHz sampling frequency. As soon as the threshold is high enough to ensure that the 
probability of having a false alarm in the analysis window is below one (i.e. for a threshold around 3), this curve is completely 
independent on the value of N used in the simulation. 
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FIG. 2. Effect of the consecutive filter output correlation on the false alarm rate. The top graph compares the normalized 
false alarm rate Tnorm for a single detector - as computed with the simulation procedure described in the core of the paper 
- with the false alarm rate rindcp computed for the same threshold, assuming that successive filter outputs are uncorrelated 
normal random variables. The real false alarm rate is at least one order of magnitude below the Gaussian estimate, on the 
whole range of threshold considered. The continuous line in the bottom plot shows, as a function of the threshold r/, the mean 
size (in bins) of the false alarm clusters which is, to a very good approximation, equal to the ratio rindep / TViorm represented by 
the black bullets. 
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Single interferometer 
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FIG. 3. ROC for a single interferometer with u) = 1 ms. Five curves corresponding to different values of the maximal SNR 

- Pmax = 5, 7.5, 10, 12.5 and 15 - are plotted. Like for all the forthcoming ROC, vortical lines give convenient conversions of 
the false alarm rate per bin, assuming a sampling frequency /samp ~ 20 klfz; from left to right, one false alarm per year, per 
week, per day, and per hour respectively. Comparing these curves clearly show that one cannot have a high detection efficiency 
associated to an high confidence level - i.e. a very small false alarm rate - by using only one detector to seek GW bursts. 
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P 

FIG. 4. Timing resolution A<;rms (in ms) of the Gaussian filter versus the SNR p of the signal as detected in the inter- 
ferometer. When p ^ 0, the noise becomes dominant: the timing of the maximum output is uniformly distributed in the 
analysis window (N = 512) and the resolution reaches the plateau N/vT^//sampiing ~ 7.4 ms. The analytical fit presented in 
[14] AtRMS = 1.45/p is valid in the intermediate range p S [6; 30]. For smaller values of the SNR, it underestimates the timing 
uncertainty which increases much faster because the noise contribution becomes more dominant. On the other hand, the fit 
overestimates AtRMS for very high values of p. 
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Virgo-LIGO loose twofold coincidences 
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FIG. 5. ROC comparing the pairs of detectors belonging to the three-interferometer network Virgo + the two 4 km-LIGO 

detectors. In this graph, the maximum SNR, pmax is set to 10. As the two LIGO interferometers have been built together in 
order to be 'aligned', this network shows bettor performance than the two other ones. 
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LIGO loose coincidences 
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FIG. 6. ROC characterizing the two LIGO detectors which form the best pair of interferometers. The five curves cover the 
range between pmax = 5 and 15. Even for the larger values of the optimal SNR, detection efficiency remains below 50% for 
manageable false alarm rates (below 1/s). 
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FIG. 7. ROC comparing all the coincidence strategies in the Virgo-LIGO network for pmax = 10. The best configuration 
requires at least two detections among three, but its efficiency remains limited for a false alarm rate small enough; in addition, 
the full coincidence appears unlikely. These two results clearly show that larger networks are required. A last point worth 
being mentioned is that for manageable false alarm rates, the detection efiiciency is better for a two interferometer network 
than for a single detector. 
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Six-interferometer loose coincidences 



73 



Q 
4= 



O 



5fold 



2fold coincidences, • 



3fold coincidences 



coincidences 



Full 



coincidence 



o 



4fold coincidences 



^ _-|2 -11 -10 -9 -8 -7 -6 -5 -4 

10 10 10 10 10 10 10 10 10 



F.A. rate per bin 



FIG. 8. ROC for the full network coincidence strategies, from twofold (top curve) to sixfold coincidences (bottom curve) 

and with pmax = 10. Detection efficiencies clearly improve by going to three to six interferometers; indeed, in this larger 
configuration, both twofold and threefold coincidences are likely, even at very low false alarm rates. 
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Virgo-LIGO and LIGO network loose coincidence comparison 
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FIG. 9. Comparison between two networks of three interferometers: Virgo and the two LIGO 4-km detectors on the one 
hand, and the three LIGO detectors (including the LIGO-Hanford 2-km) on the other hand. In order to simplify the labels of 
the two plots, the interferometer names are shortened: Virgo (V), LIGO Hanford 4-km (H4km) and 2-km (H2km) and LIGO 
Livingston 4-km (L4km). The top graph presents ROC corresponding to coincidences between pairs of detectors: from top to 
bottom, the two LIGO 4-km detectors, Virgo associated with each of the two LIGO 4-km interferometers, the two Hanford 
detectors (4 km and 2 km) and finally Livingston 4-km with Hanford 2-km. The bottom graph compares the two possible 
strategies involving all the detectors of these networks: twofold coincidences (at least two detections among three) and full 
coincidences. As the two LIGO Hanford detectors have identical locations, their coincidence false alarm rate is lower than 
for any other pair of detectors. Yet, this docs not compensate the difference in sensitivity between them which limits their 
detection efficiency: the full LIGO network is less efficient than the Virgo-LIGO detector. This clearly shows the importance 
of the final interferometer sensitivity. 
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Tight coincidences comparison 
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FIG. 10. Comparison of ROC (pmax = 10) corresponding to various tight coincidence strategies: LIGO coincidences, twofold 
and threefold detections in the Virgo-LIGO network and in the full network of six interferometers. 
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FIG. 11. Sky maps and distributions of | cosOab] as a function of the source sky coordinates {a, 5) for two different network 
configurations: (top) Virgo and the two LIGO interferometers ; (bottom) the full set of six antennas. 
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Virgo-LIGO coherent data analysis 
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FIG. 12. ROC for the coherent search of GW burst signals in the three-interferometer network Virgo-LIGO, assuming the 
source sky location to be known. The curves have been computed for four different values of the optimal SNR: pmax = 1, 5, 
7.5 and 10 respectively. 
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Six-interferometer coherent data analysis 
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FIG. 13. ROC for the coherent search of GW burst signals in the full network of interferometers including the 6 currently 

existing projects in the world. For pmax = 10, the detection efficiency is higher than 95% in the whole range of r considered 
and e > 80% for the case pmax = 7.5, assuming the source location to be known. 
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Virgo-LIGO network 
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FIG. 14. Timing resolution of the coherent analysis showing the evolution of the timing error AtRMS (in ms) as a function 

of the square root of the coherent statistics (Aroducod)^^'^. The two curves for the Virgo-LIGO network and for the full set of 
interferometers - are identical, as expected from the network- independent statistics Areduced. 
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